Kondo effect of non-magnetic impurities and the co-existing charge order in the 

cuprate superconductors 

Kwon Park 

Department of Physics, Yale University, P.O. Box 208120, New Haven, CT 06520-8120 

(February 1, 2008) 

We present a theory of Kondo effect caused by an induced magnetic moment near non-magnetic 
impurities such as Zn and Li in the cuprate superconductors. Based on the co-existence of charge 
order and superconductivity, a natural description of the induced moment and the resulting Kondo 
effect is obtained in the framework of bond-operator theory of microscopic t-J-V Hamiltonian. The 
local density of state near impurities is computed in a self-consistent Bogoliubov-de Gennes theory 
which shows a low-energy peak in the middle of superconducting gap. Our theory also suggests that 
the charge order can be enhanced near impuries. 
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Impurities in the cuprate superconductors have at- 
tracted much interest because of their role as an effective 
tool in probing the local electronic and magnetic prop- 
erties of high Tc superconductors. Especially remarkable 
have been recent experiments in scanning tunneling mi- 
croscopy (STM) of surfaces of the cuprate super- 
conductor Bi2Sr2CaCu208+5 (BSCCO) which showed a 
low-bias peak in the differential tunneling conductance in 
the vicinity of Zn impurities Since the Zn++ ion has 
spin 5 = 0, it seems natural to view each Zn site as just 
a potential scatterer, and therefore to interpret the peak 
as a quasibound state However, there are some ex- 
perimental features which cannot be naturally explained 
by this theory Q| . 

More importantly, a series of NMR experiments [^|J^] 
have shown clear evidences that there is an induced mag- 
netic moment with spin S = 1/2 near Zn impurities. 
Therefore, it is imperative for a consistent theory to ad- 
dress the problem of impurity in the framework of Kondo 
effect. Though there have been several, previous works 
based on the Kondo physics their approach was 

based on the ad hoc assumptions on the existence of in- 
duced spin moment near Zn impurities. Especially, the 
location of induced spin moment and the range of Kondo 
interaction were chosen in ad hoc manner. In this article, 
we would like to provide a self-consistent theory of Kondo 
effect based on the co-existence of charge order and su- 
perconductivity. We take this assumption to be natural 
because one of the important consequences of charge or- 
der is the induced S = 1/2 moment near non- magnetic 
impurities as schematically shown in Fig.|^. 

Much more direct evidences for the relevance of charge 
order in the high Tc superconductors have been recently 
obtained in several remarkable STM experiments, one of 
which was performed by Hoffman et al. |^ where the 
charge order with a period of four lattice spacing was 
observed near vortex cores in BSCCO. As predicted in 
Ref. [D, the superconductivity is locally suppressed near 
the cores of vortices, and the static charge order ap- 



pears in such regions. Another STM measurement in 
optimally doped BSCCO was performed by Howald et 
al. [|lO|, showing clear evidences for the static charge 
order even in zero magnetic field. In addition to the 
STM measurements, inelastic neutron scattering tech- 
niques have been also used to provide important evi- 
dences for the charge order in La2-xSr2.Cu04 (LSCO) 
and YBa2Cu306+:. (YBCO) |n|,|l|. 

Note our viewpoint that, though the charge order may 
not be static in all regions of superconducting phase, the 
time scale of charge order fluctuation is long enough to 
ensure the superconducting correlation and the Kondo 
spin physics. Therefore, in the rest of discussions, we will 
assume the static charge order with period of two lattice 
spacing (the spin-Peierls order), which, we believe, cap- 
tures the essential physics in spite of its simplicity. Also, 
remember that the charge orders with different period 
are most likely to coexist because, as shown in Ref. 
the periodicity is rather sensitive to various parameters 
including the doping concentration. 

Recently, a natural theoretical framework was devel- 
oped to describe the co-existing phase of charge density 
wave (CDW) order and superconductivity, based on the 
bond-operator representation It was shown that the 
saddle-point approximation of this bond-operator the- 
ory consistently interpolates the phase of Mott insu- 
lator at low dopings and that of superconductor with 
nodal fermions at moderate dopings. Therefore, it would 
be very interesting to explore the question whether the 
bond-operator theory gives rise to the low-energy peak 
of local density of state in the vicinity of non-magnetic 
impurities without ad hoc assumptions. It would be also 
interesting to study the effect of impurity on the charge 
order. 

Now we begin our bond-operator theory by setting up 
the exact mapping between the bond operators and the 
usual electron creation operators [||. Let c\^ and c^^i 
(a =Tii) be the electron creation operators on the two 
sites of a pair. When we project out all states with two 



electrons at the same site, the electronic Hilbert space 
for a pair of sites is composed of nine states which can 
be expressed in terms of the "bond particle" creation 
operators defined by: 
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where |0) is the electron vacuum and \v) is an imagi- 
nary vacuum void of any bond particles. Remember that 
ct"j (a — x,y,z) are the Pauli matrices, and Eat is the 
second-rank antisymmetric tensor with £|| = +1. The 
opertators s,d,ta all obey the canonical boson commu- 
tation relations, while the /iia,^2a obey the canonical 
fermion relations. Since the total Hilbert space of these 
five bosons and four fermions is much larger than that 
of the physical nine states, we must impose the following 
constraint on the bond particle Hilbert space. 



s + tita + h\^hia + hl^h2a + dU = 1. (6) 



In the subspace contrained by Eq. (|6|) we can write the 
exact expressions for electron operators in terms of bond 
operators: 



c\a = h\ad + -^eafoS^/lafc - -^eac(Tcbtih2b: (7) 

4q = ^Aad + -^eabS^hu + -^EacfJcbtihib, (8) 

Sla ^{s'^ta + tls - eafi^t^pt^) + ^CTab^lahlb, (9) 

= -7;is''ta + + eap-At^) + \Kbh\ah2b, (10) 



where tap-y is the third-rank antisymmetric tensor with 
^xyz = +1, and, as usual, the electon spin operator is 
defined as Sa = l/2c\a'^f^Cb- 

Equiped with the bond-operator representation of elec- 
tron operators, we now apply the Bogoliubov-de Gennes 
(BdG) theory to solve the t-J-V model in finite systems 
containing a single impurity with (infinitely) strong re- 
pulsion, J7, which is defined by: 

H = ^(Cja'^JQ + ^a^ia) + ^ ^ SiaSja 
(id) (id) 

+ ^ H 4aCmc]f,Cj6 '^'^^Cm + Ucl^d^a (H) 

i^d) i 

where ia denotes the position of impurity, and fi is the 
chemical potential. Also, it is implicitly assumed that all 
states with two electrons on any site have been projected 



out. Written in terms of bond operators, the Hamilto- 
nian is physically meaningful only when the contraint 
condition in Eq.(^ is satisfied simultaneously. In this 
article, as usual in the mean field theory, only the aver- 
age of this constraint will be satisfied via the Lagrange 
multiplier method. Note, however, that the average is 
taken over quantum fiuctuations, not over space, so that 
the constraint will be satisfied individually at each site. 

While general techniques of the bond-operator method 
can be found in Ref. (ij] in detail, several conceptial 
and technical points are worth mentioning, especially re- 
lated to the impurity problem. (1) The nearst-neighbor 
Coulomb repulsion, V, is not only physically reasonable, 
but also gives rise to an important consequence. With- 
out V, the pairing of holes primarily occurs through the 
condensation of d-bosons, which results in a very short- 
range s-wave-like pairing. It is only when the d-boson 
condensation is suppressed by a large V that the super- 
conducting state develops c?-wave-like pairing, and nodal 
fermions emerge. Remember that we do not make any 
assumptions either on the emergence of superconducting 
state or the symmetry of pairing. They are obtained as 
a natural consequence of our saddle-point bond-operator 
theory of t-J-V model at moderate dopings. 

(2) The chemical potential, fi, should be determined 
by fixing the average hole concentration, x, in the region 
far away from the impurity. In our finite system, we first 
obtain the value of /x without introducing the impurity, 
and then use it for the case of impurity by assuming that 
the system size is large enough so that a single impurity 
does not change the chemical potential. 

(3) The impurity potential, U, is taken to be infinitely 
repulsive so that electrons are completely depleted from 
the Zn site. In the bond-operator formalism, one of the 
/i-fermions, say hia{a =t,i), is pinned at the unpaired 
site near the Zn impurity. (See Fig.^) Of course, it is 
possible without charge order in the doped antiferromag- 
nets that the electron escapes from the unpaired site all 
together, leaving the vacancy instead of lone spin mo- 
ment. It is assumed, however, that the empty state is 
energetically unfavorable because of the similar reason 
why d-boson condensation is suppressed for large V . 

(4) Finally, unlike the previous work in Ref. Q, there 
are now two sets of BdG equations both for the bosons 
{ta) and the fermions {hia and /i2a)- Consequently, the 
number of self-consistency conditions for the normal and 
anomalous exchange energies is greatly increased, which 
limits the system size significantly. The number of self- 
consistence equations is 11 x with N'^ being the 
number of sites. Also, note that the diagonalization 
of boson BdG matrix is not as straightforward as that 
of fermion BdG matrix. It is performed by using the 
method in Ref. [H). 

The saddle-point theory is not yet complete alone with 
the aforementioned self-consistent BdG equations, but it 
still needs to determine the condensation density of s- 



boson, (s) , by minimizing the ground state energy. Sim- 
ilar to the chemical potential, the Lagrange multiplier as- 
sociated with the constraint is first computed as a func- 
tion of (s)^ for the situation when there is no impurity, 
by satisfying the constraint in the mean field level. Then, 
(s)^ is fixed to be the value minimizing the ground state 
energy, at which point the Lagrange multiplier is also 
fixed. Fig.^ shows the ground state energy as a function 
of s-boson condensation density which also can serve as 
an order parameter for the spin-Peierls order. After the 
impurity is introduced, the position-dependent conden- 
sation density, (s^)^, is determined in turn by satisfying 
the constraint at each site with the Lagrange multiplier 
determined previously. In this way, it is guaranteed (at 
least in the mean field level) that the ground state en- 
ergy is minimized and also the constraint is satisfied at 
the same time. 

Now let us turn to the physical observables. As men- 
tioned in the introduction, the local density of state 
(LDOS) can be obtained directly in the STM experiments 
through the measurement of differential tunneling con- 
ductance, which show a low-energy peak in LDOS near 
non-magnectic impurities while the superconducting co- 
herence peak is greatly suppressed Q . Theoretically, the 
local density of state can be computed in terms of the 
spectral function given by: 

p(r,w) =Im (r,w + i<5), (12) 

where the retarded electronic Green function, G"''^*(r, cj + 
iS), is obtained by using the usual analytic continu- 
ation from the time-ordered electronic Green function 
G(r,r) - (r,Ca(r,T)ct(r,0)) (a =T,i). In the bond- 
operator formalism, the spectral function of electron on 
the site 1 of dimer located at is related to that of /12- 
fermion via: 

pi{r,,Lo) - -^(s,f Im G-^f (r„ - z6), (13) 

where Gh2{iCi,T) = (r^/i2a(ri, T)/iJ„(ri, 0)). The elec- 
tronic spectral function at site 2 is related to that of 
/ii-fermion in similar way. In Eg. (p^) the approximation 
is made when the contribution from magnons (ta) is ig- 
nored because they are high-energy modes. 

To compare the theoretical LDOS with the differential 
tunneling conductance measured in STM experiments, 
we need to take into account some key properties of re- 
alistic surface structure of BSCCO. Between the super- 
conducting Cu02 layer and the STM tip, there is always 
the BiO layer formed in such a way that each Bi atom 
is located directly above each Cu or Zn atom. So it is 
reasonable to assume that the Bi atom will block tun- 
neling currents from reaching the Cu/Zn atom directly 
below the STM tip, and so the STM measures the LDOS 
contributed by the four nearest-neighboring sites instead 
of the single site directly below the tip. Under 



this assumption, the LDOS measured in STM experi- 
ments, {p{vi,uj)), may be identified with an average of 
LDOS from the nearest-neighboring sites: {p{ri,uj)) (x 
p{ri + x,uj) + p{vi ~ x,uj) + piji + y,uj) + p{vi - y,uj). 
Fig.^ shows {p{Yi, E)) near and far away from the Zn im- 
purity as a function of energy E/ J, where the LDOS near 
impurity is defined as {pivi^uj)) with r; indicating the 
position of the Zn impurity. Compared to the LDOS far 
away from the impurity, the low-energy peak (denoted 
by the arrow in graph) develops near the Zn impurity, 
while the coherence peak is somewhat reduced indicat- 
ing the suppression of superconductivity, which is in a 
reasonably good agreement with STM measurements. 

Another physical observable in which we are interested 
is the order parameter of charge order, or the spin-Peierls 
order parameter, (s^) . One of the main questions that 
we would like to answer in this article is how much the 
charge order is affected by the impurity. It is completely 
conceivable that the spin-Peierls order is locally reduced 
near impurity at the same time when the superconduc- 
tivity is also suppressed. In the extreme case the spin- 
Peierls order can be completely destroyed while the 5=1 
excitons start to condense locally so that the magnetic 
order emerges near the impurity: in the bond-operator 
formalism, {ta{vi)) 7^ with indicating the neighbor- 
ing sites of impurity. If so, it may imply the emergence of 
the spin density wave (SDW) near impurities. However, 
the {si) plotted as a function of the distance from im- 
purity in Fig.^ shows the opposite behavior that overall 
the spin-Peierls order is not affected much by the pres- 
ence of impurity after one lattice spacing. In fact, the 
spin-Peierls order is somewhat enhanced at the nearest- 
neighboring Cu site of the Zn impurity, which may sug- 
gest an interesting prediction that the static charge can 
be observed near the Zn impurities by using the STM 
techniques. Finally, note that the bottom panel of Fig.^ 
shows the equal-time spin-spin correlation within dimer, 
^ii^2i = "I'^i + i4*i' which suggests basically the sim- 
ilar physics by showing that the effect of the spin S = 1 
excitons is small. 

In conclusion, we have applied the bond-operator 
method to address the problem of non-magnetic impu- 
rity such as Zn/Li in the high- Tc cuprate superconduc- 
tors. It is shown that, without ad hoc assumptions, the 
low-energy peak of local density of state in the vicinity 
of impurity can be computed in the framework of bond- 
operator formalism. Also, it is predicted that the charge- 
density-wave order can be enhanced near Zn impurities. 

The author is indebted to Subir Sachdev for very use- 
ful discussions and encouragement. This research was 
supported in part by US NSF Grant DMR 0098226. 
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FIG. 1. Schematic diagram of the ground state of doped 
antiferromagnet with Zn impurity, which is assumed to have 
bond-centered charge order with period of two lattice spac- 
ing (columnal spin-Peierls order). Ellipses denote the valence 
bonds formed by two neighboring spins: -^(| t) i) ^ I i) T))- 
The arrow near Zn impurity indicates the induced spin mo- 
ment. Note that a similar picture of moment formation ap- 
plies to other charge order states. 
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FIG. 2. Ground state energy as a function of the conden- 
sation density of s-boson for the sqaure lattice system of 8 x 8 
sites without the impurity. Here t/ J = 1.5 and the hole con- 
centration X = 0.3. 
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FIG. 3. Local density of state (LDOS) near the impurity 
(solid histograms) and far away from the impurity (dashed 
histograms) as a function of energy in the finite square-lattice 
system of 12 x 12 sites with t/ J = 1.5 and x = 0.3. The dashed 
line is just a guide to eye for the LDOS far away from impurity, 
whose shape is obtained from the analytic computation for 
the uniform system, and is fitted roughly to follow the LDOS 
computed in finite system. Remember that the LDOS is not 
normalized due to the discrete nature of energy spectrum in 
finite system. It is important to note that the low-energy 
peak (denoted by arrow) develops in LDOS near the impurity, 
while the coherence peak is suppressed, which is in a good 
agreement with experiments. 
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FIG. 4. Condensation density of s-boson, (si)^ (top 
panel) , and the equal-time spin-spin correlation within dimer, 
(SfiSfi), (bottom panel) as a function of distance from the 
impurity in the finite square-lattice system of 12 x 12 sites with 
t/J = 1.5 and x = 0.3. Remember that in the bond-operator 
formalism SiiS2i = — fs? -I- jtlti whose expectation value 
equals to —3/4 for the pure singlet valence bond, and 1/4 for 
the pure triplet magnon. Of course, both the condensation 
density and the correlation are exactly zero at the impurity 
site since the Zn impurity is assumed to avoid electrons com- 
pletely. Circles (crosses) in the graph are associated with the 
direction parallel (perpendicular) to the columnal spin-Peierls 
order. It is interesting to observe that the spin-Peierls order, 
or the charge density wave, is somewhat enhanced near the 
impurity while overall it is not affected much by the presence 
of impurity after one lattice spacing. 



